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ABSTRACT 

■ We provide scaling relations and fitting formulae for adiabatic cold dark matter cosmologies 
that account for all baryon effects in the matter transfer function to better than 10% in the large- 
scale structure regime. They are based upon a physically well-motivated separation of the effects 
of acoustic oscillations, Compton drag, velocity overshoot, baryon infall, adiabatic damping, 
Silk damping, and cold-dark-matter growth suppression. We also find a simpler, more accurate, 
and better motivated form for the zero baryon transfer function than previous works. These 

i descriptions are employed to quantify the amplitude and location of baryonic features in linear 

theory. While baryonic oscillations are prominent if the baryon fraction fib/fio J> £loh 2 + 0.2, 
the main effect in more conventional cosmologies is a sharp suppression in the transfer function 
below the sound horizon. We provide a simple but accurate description of this effect and stress 
£NJ | that it is not well approximated by a change in the shape parameter T. 
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■ Subject headings: cosmology: theory - dark matter - large-scale structure of the universe 
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jH . 1. Introduction 

A key success of the cold dark matter (CDM) paradigm is the ability of linear perturbation theory 
in the early universe to explain the power spectra observed in the cosmic microwave background (CMB) 
and galaxy surveys. On the largest scales, COBE finds a trend of power with scale rather close to the 



theoretically motivated scale- invariant spectrum (Bennett et al. 1996). On scales between 10 and 200 Mpc, 



however, galaxy surveys (e.g. Baugh & Efstathiou 1993; Btrauss fc Willie^ 1995) find a much different 
trend, in which the power increases with scale. Merely by including the effects of the transition between a 
radiation-dominated and matter-dominated universe, the CDM cosmology roughly explains both the relative 
normalization and the differing spectral indices of these two regimes. 

Although the presence of cold dark matter does play a leading role in determining the matter power 
spectrum, the inclusion of baryons can lead to significant alterations. Indeed baryonic features in the power 
spectrum are a fundamental prediction of the gravitational instability paradigm, and their discovery would 
represent a strong consistency test for the cosmological model. Such features are the direct result of small 
density fluctuations in the early universe prior to recombination. At those times, the baryons are tightly 
coupled with the photons and share in the same pressure-induced oscillations that lead to acoustic peaks in 
the CMB. This not only leads to intermediate-scale oscillations in the power spectrum but also produces an 
overall suppression of power on small and intermediate scales. 
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While the low baryon fraction (~ 5%) in the standard CDM model may have justified the neglect of 
baryonic effects on the power spectrum in the past, recent observations favor higher baryon fractions. X-ray 
observations of clusters of galaxies yield baryon fractions of 10-30% (White et al. 1993; David et al. 1995). 



Recent measurements of high-redshift deuterium abundances ( Tytlcr et al 



1996; but see 



Rugers fc Hogan 



1996) and new theoretical interpretations of the Lyman-a forest ( |Weinberg et al. 1997, and references therein) 
suggest a value of the baryon density f2f, greater than the fiducial nucleosynthesis value of 0.0125/i~ 2 (Walker 
1991). Meanwhile, observations of large-scale structure (c.f. Bahcall| 1996), higher Hubble constants 



et al 



(c.f. Frccdman 1996), and high-redshift galaxies (e.g. Spinrad et al. 1996) and clusters ( Bahcall et al. 1997) 
favor a universe with Qq < 1. Such baryon fractions lead to modifications of the pure CDM transfer function 
that are within observational sensitivities ( Tegmark 1997; |Goldbcrg fc Strauss) 1997). 



Numerical codes to solve the multi-species Boltzmann equations (e.g. Bond fc Efstathiou| 1984; Holtz 



manr 1989; |Hu et al.| 1995; |5eljak fc Zaldarriaga] 1996) now agree to 1% accuracy and run in a few minutes on 
today's workstations. While these codes (e.g. the publically available CMBfast) should be used for applica- 
tions demanding high accuracy, analytic descriptions are useful for understanding how the different physical 
effects give rise to the behavior seen in the transfer function. Such descriptions better isolate the unique and 
robust observational signatures of physical processes in the early universe and quantify their scalings with 
cosmological parameters. They also probe possible parameter degeneracies and suggest possible consistency 
tests with related effects in the CMB. 

To this end, we develop here a fitting formula for the matter transfer function of the general CDM 



plus baryon universe [c.f. eqs. (16)-(|2j) and §2]. The formula is composed of a number of well-motivated 
ingredients, whose behaviors we discuss in detail. We achieve fractional accuracies of 10% in fully-baryonic 
models and ~ 5% in partial baryon models. Included here is a quantification of the fundamental scales 
including the acoustic and Silk damping scales which are related but not equal to the equivalent scales in 
the CMB. 

We then use this form to produce quantitative assessments of the amplitude and location of baryonic 
oscillations as well as the alteration to the intermediate-scale shape and small-scale normalization of the 
transfer function. Previous assessments of the latter effects (see e.g. Peacock fc Dodds 1994; Bugiyama 



1995), though sufficient for the current generation of experiments at low baryon fraction, are not accurate 
enough for future high-precision tests expected of the Sloan Digital Sky Survey ( Gunn fc Weinberg] 1995) and 
2-Degree Field survey ( Taylor! 1995). We present an approximate form that neglects the acoustic oscillations 
but accurately represents the suppression of power on intermediate scales [c.f. eqs. ( p6| ) and (28)-(|3l|)]. To 
this end, we also present a simpler and more accurate formula [eq. (p9|)] for the zero-baryon transfer function 
(e.g. Bardeen et al. 1986). 



In §2, we lay the groundwork for the subsequent discussions by presenting a summary of the physical 
scales that enter linear perturbation theory and the exact small-scale solution (Hu & Sugiyama 1996, here- 
after HS96) that we use to anchor our fitting formula. In §3, we state the fitting formula and discuss its 
performance. In §4, we describe the phenomenology revealed by the formula and present simple scalings to 
characterize the baryon oscillations and shape alteration. We conclude in §5. In Appendix A, we give a 
short guide to help the reader locate relevant formulae from the paper and turn them into COBE-normalizcd 
power spectra. A list of symbols used in this paper is given in Table |. 
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Effective shape 
Temperature of CMB 
r suppression 
Baryon suppression 
CDM suppression 
Baryon envelope shift 
Shift in CDM log- 
Node shift parameter 
Photon-baryon ratio 
R at Zd 

Transfer function (TF) 

Baryon sector TF 

CDM sector TF 

Zero baryon TF 

Pressureless TF 

Wavenumber 

Equality wavenumber 

Silk wavenumber 

k scaled with fc cq 

Sound horizon at drag epoch 

Effective sound horizon 

Redshift of drag epoch 

Redshift of equality 




Table 1: A list of symbols used in this paper. 



2. Physical Effects 

To motivate and explain the transfer function formulae in §[|, we begin with a review of some of the 
basic results of linear perturbation theory, starting with a summary of the physical scales that enter the 



theory. We then describe the exact small-scale solutions found by HS96, as they play a central role in the 



development of the fitting formulae and the explanation of baryon phenomenology. 

The particular physical properties of the constituents of the universe, in particular their equations of 
state and interactions, can alter the predictions of perturbation theory. Causality, however, precludes effects 
at arbitrarily large scales. It is therefore usual to measure the resulting perturbations by comparing them 
to the amplitude they would have had were causal physics neglected. The result is the transfer function, 
defined asp] 

T(b\- 6 ( k ' z = °) a(0,* = oo) 
1 >- S(k,z = oo) (5(0, 2 = 0) ' 1 ' 

where 5(k,z) is the density perturbation for wavenumber k and redshift z. By construction, T — > 1 as 



2 in the synchronous or comoving gauge; gauge issues of are no practical concern for subhorizon scales. 
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k — ► 0. The power spectrum (|<5(k)| 2 ) is proportional to the square of the transfer function multiplied by 
the initial power spectrum, most often taken to be proportional to a power law k n with n ~ 1. Strictly 
speaking, each species of particle has a separate transfer function; however, after recombination the baryons 
are essentially pressureless and quickly catch up with the cold dark matter perturbations, leaving both with 
the same transfer function. It is this transfer function that we study in this paper. 

We consider cosmologies in which the universe is primarily composed of photons, baryons (and their 
accompanying electrons), massless neutrinos, and cold dark matter (CDM). Relative to the critical density, 
the densities today of the baryons and CDM are fib and fl c respectively. The total matter density is then 
flo = fib + fl c . The CMB temperature Tqmb is written as 2.702.7 K; the best observations from the COBE 



FIRAS instrument are 2.728 ± 0.004 K (Fixsen et al. 1996; 95% confidence interval). We assume that the 



massless neutrinos contribute an energy density corresponding to three species at (4/11) 1 / 3 the temperature 
of the photons. We use a Hubble constant H and define h = iJ /(100kms _1 Mpc -1 ). It is important to 
note that for z ^> Qq 1 the dynamics of the fluctuations are determined solely by the matter-radiation ratio 
5!o/i 2 02~7, the baryon-photon ratio flbh 2 ^^ an d the CMB temperature Xcmb- We fix the last of these to 
be the COBE value and do not include variations in 02.7 in our fits. Since all effects in the transfer function 
are set at those early times, the resulting description should depend only on floh 2 and flb/flo- The existence 
today of a non-zero cosmological constant or curvature is insignificant. 



2.1. Length and Time Scales 

The physics governing the evolution of perturbations in CDM-baryon universes involves three distinct 
length scales: the horizon size at matter-radiation equality, the sound horizon at the time of recombination, 
and the Silk damping length at recombination. 

In the usual cosmological paradigm, non-relativistic particles (baryons, electrons, and CDM) dominate 
relativistic particles (photons and massless neutrinos) in density today. However, because the density of 
these two classes of particles scale differently in time, at an earlier time, the reverse situation held. The 
transition from a radiation-dominated universe to a matter-dominated one occurs roughly at 

z cq = 2.50 x io A n h 2 e^, (2) 

the redshift where the two classes have equal density. As density perturbations behave differently in a 
radiation-dominated universe versus a matter-dominated one due to pressure support, the scale of the particle 
horizon at the equality epoch z eq , 

k cq = (2n a H%z cq ) 1/2 = 7.46 x 10" 2 fi ^ 2 e^7 Mpc" 1 , (3) 

is imprinted on the matter transfer function; in particular, perturbations on smaller scales are suppressed in 
amplitude in comparison to those on large scales. If the universe consisted only of non-interacting matter 
and radiation, the matter transfer function would depend on the ratio (k/k eq ) alone. 

Complications arise due to interactions between the species. Prior to the recombination of baryons 
and electrons, the large density of free electrons couples the baryons to the photons through Coulomb and 
Compton interactions so that the three species move together as a single fluid. This continues until, in the 
process of recombination, the rate of Compton scattering between photons and electrons becomes too low, 
freeing the baryons from the photons. We thus define the drag epoch za as the time at which the baryons are 
released from the Compton drag of the photons in terms of a weighted integral over the Thomson scattering 
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(a) Equality-Sound Horizon (b) Sound Horizon-Silk Scale 
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Fig. 1. — Comparison of the physical scales as functions of Vt$h 2 and the baryon fraction Q(,/fio- (a) The 
equality scale vs. the sound horizon: k cq s/ir (unlabeled contours at 0.1 increments), (b) The sound horizon 
vs. the Silk scale: /csiiks/7r (unlabeled contours 2 and 3). The factors of it have been included to facilitate 
comparison with the acoustic scale. 



rate (see 



HS96] , Eqn. C8, E2). A fit to the numerical recombination results is 

(n h 2 ) - 251 



z d = 1291- 



[i + fei(fW) 



T + 0.659(tt /i 2 ) - 828 ' 
b x = 0.313(r2o/i 2 )- - 419 [l + 0.607(f7 /i 2 ) a674 ] 
b 2 = 0.238(ft /i 2 ) a223 , 



(4) 



where we have reduced z c i by a factor of 0.96 from HS96 on phenomenological grounds. For Qi,h 2 < 0.03, 
this epoch follows last scattering of the photons. 

Prior to z,i, small-scale perturbations in the photon-baryon fluid propagate as acoustic waves. The 
sound speed is c s = 1/ y/3(l + R) (in units where the speed of light is unity) , where R is the ratio of the 
baryon to photon momentum density, 



i?= 3p 6 /4p 7 = 31.51We.^(z/l(r) 



3\-l 



(5) 



We define the sound horizon at the drag epoch as the comoving distance a wave can travel prior to redshift 

Zd, 



t(z d ) 2 ry 

c s {l + z)dt= — J — In 



VTTR^+^R d + R. 



cq 



eq 



CC| 



1 



(6) 



eq 
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where Rd = R(zd) and R cq = R(z cq ) are the values of R at the drag epoch and epoch of matter-radiation 
equality, respectively. The sound horizon at the drag epoch (hereafter simply the sound horizon) is larger 
than the equality horizon (~ l/fc eq ) in high-Jig models but smaller than it in low-il models; it also decreases 
strongly with increasing baryon fraction if Rj :> 1 (see Fig. |l|). 

On small scales, the coupling between the baryons and the photons is not perfect, such that the two 



species are able to diffuse past one another ( Silk 1968). The Silk damping scale is well fit by the approximation 



fcsuk = l.6(n b h 2 )°- 52 (n h 2 ) - 73 [1 + (10.4^/j 2 )- - 95 ] Mpc- 1 , (7) 



which represents a ±20% phcnomcnological correction from the value given in HS96 . The Silk scale is 
generally a smaller length scale than either s or l/k eq . Note that the difference between the drag and last 
scattering epochs implies that for Sl^ti 2 <; 0.03 the sound and Silk scales in the transfer function are larger 
than those in the CMB. 

We show a comparison of these scales as a function of cosmological parameters in Figure |[ 



2.2. Small-scale Solutions 

In the small-scale limit, one can solve the linear perturbation equations analytically in the approximation 



that baryons provide no gravitational source to the CDM ( HS96 ) . This approximation is appropriate below 
the sound horizon since baryon perturbations are pressure supported. As we will use this solution in order 
to anchor the small-scale end of our fitting formulae, we describe the solutions further. 

The transfer function is written as a sum of the baryon and cold dark matter contributions at the drag 
epoch 

T(k) = ^T b (k) + ^T c (k). (8) 

The CDM transfer function can be solved exactly in terms of hypergeometric functions that are more 
conveniently approximated by the following form: 



«c^^ (9) 



In 1.8p c q 
U.2q 2 

(fc/Mpc-^e 2 .^ 2 )- 1 ^— k — , (10) 



13.41fc oq 



where a c and (3 C are fit by 



_ -Jib/Oo -(fijj/fio) 3 

ai = (46.9^ 2 )°- 670 [l + (32.m / l 2 )- - 532 ], 

a 2 = (12.0r!o^ 2 ) a424 [l + (45.0fi / l 2 )- - 582 ], 

0- 1 = l+M(^cM)) b2 -i], (12) 

h = 0.944[1 + (458fto/i 2 r a708 r\ 

b 2 = (0.395fW)-° 0266 . 

As fib/f^o — * 0, a c ,Pc — > 1. Equation (||) shows the familiar ln(fc)/fc 2 dependence of the small-scale CDM 

transfer function. This occurs because outside the horizon, density perturbations grow as k 2 due the product 




of potential and velocity gradients which drive the growth; inside the horizon in the radiation dominated 
epoch the growth is logarithmic. The main effect of the baryons comes from the suppression in growth rates 
between equality and the drag epoch. As floh 2 increases the time between the two epochs increases; thus the 
maximum suppression due to the baryons occurs in the highest Q h 2 models. A plot of a c Q c /Slo is shown 
in Figure |^a. 

In the small-scale limit, the baryons are trapped in acoustic oscillations until recombination permits 
them to slip past the photons. While the density perturbation of this oscillation contributes to the transfer 
function, the corresponding velocity perturbation actually dominates in the small-scale limit. When the 
oscillation is released at the drag epoch, the baryons move kinematically according to their velocity and 



generate a new density perturbation (Sunyaev & Zel'dovich 1970; Press & Vishniac 1980). This so-called 



velocity overshoot means that the transfer function for ks ^> 1 follows 

Tb -> ab S ^hflv(k). (13) 
ks 

Here T>(k) represents the effects of Silk damping, which occurs due to combination of diffusion of the photons 
with respect to the baryons and Compton drag moving baryons from overdensities to underdensities and 
hence destroying the perturbation. That the dependence is sin(fcs) rather than cos(fcs) is the result of the 
dominance of the velocity term rather than the density term. A detailed treatment allows one to calculate 

Olb- 

a b = 2mk cq s(l + R d )- 3 / i G(\±^) , (14) 
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G{y) 



-6yiT^+(2 + 3y) In 



( VT+y + 1 



(15) 



The factor (1 + Rd) 3 ^ 4 comes from the damping of oscillations resulting from the adiabatic decrease in 
the sound speed (|Peebles fe~Yr] 1970; [HS96| , eq. [A17]). The factor involving G{y) (oc y~ 1/2 for y > 1) 
comes from the product of the growth suppression between equality and the drag epoch (oc and the 
time available before the velocities creating the perturbation decay due to the free expansion of the universe 
(oc y 1 / 2 )- A plot of afoOb/fio is shown in Figure ||b. 

That the phase of the oscillations is ks ( |Pceblcs fc Yu 1970) can be seen simply from integrating the 
dispersion relation uj = kc s . The change in phase for an acoustic wave is 8<j) = [k(l + z)]c s St. Integrating 
this from t « to the drag epoch (where the baryons are released and the oscillations freeze-out) yields ks, 
owing to the definition in equation (j^). A technical complication occurs for k ;> fcsiik- The presence of strong 
damping slightly raises the redshift at which the oscillations freeze-out, making s a few percent smaller (see 



HS96 Fig. 2). We have neglected this effect because it occurs at sufficiently small scales that the resulting 



phase shift is unobservable in practice, but one can see the deviations when comparing to numerical results 
(c.f. Fig. |. 



3. Fitting Formulae 

As we have seen in §0, analytic solutions exist for the transfer function at both large and small scales. 
The transition between these extremes is defined by two scales, the horizon at matter-radiation equality 
and the sound horizon at the end of the drag epoch. The former sets the dynamics of the expansion and 
perturbation growth; the latter sets the scale at which pressure support becomes important for the baryons. 
Because the range of scales accessible by the study of structure formation falls within this transition regime, 
it is important to understand the full transfer function in detail. To that end, we present in this section a 
fitting formula that approximates the full transfer function on all scales. 

We write the transfer function as the sum of two pieces, 

T(fc) = ^T b (fc) + ^T c (fc), (16) 

whose origins lie in the evolution before the drag epoch of the baryons and cold dark matter respectively. This 
separation is physically reasonable as before the drag epoch the two species were dynamically independent 
and after the drag epoch their fluctuations are weighted by the fractional density they contribute. This 
automatically includes in T c the effects of baryonic infall into CDM potential wells. Note however that Tf, 
and T c are themselves not true transfer functions as they do not reflect the density perturbations of the 
relevant species today. Rather, it is their density-weighted average T(k) that is the transfer function for 
both the baryons and the CDM. 



3.1. Cold Dark Matter 



The transfer function for cosmologies in which non-interacting cold dark matter dominates over baryons 
have been studied by many authors and accurate fitting formulae already exist in this limit (e.g. Bond & 
Efstathio"u| 1984; Bardeen et al. 1986; but see improvements in §4.2) However the effect of baryons, though 
long known from numerical calculations (e.g. Peebles & Yu 1970; Holtzmann 1989), have in the past either 
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been included in fitting formulae in an ad hoc manner (see e.g. Peacock & Dodds 1994; Sugiyanra] 1995) or 
only in the small scale limit ( HS96 ). 

In the presence of baryons, the growth of CDM perturbations is suppressed on scales below the sound 
horizon. The change to the asymptotic form can be calculated and has been shown in equations ([)])-( 12). 
We introduce this suppression by interpolating between two solutions near the sound horizon: 



T c (k) 
f 



fT (k,l,p c ) + (l-f)T {k,a c ,P c ) 
1 



1 + (ks/5A) 4 



with 



T (k,a c ,(3 c ) = 



c = ™ + 

a c 



ln(e + 1.8f3 c g) 
\n(e + 1.8P c q) + Cq 2 ' 
386 



1 + 69.9g 108 ' 

q, a c , and /3 C have been given in equations (|l0|), (|TT|) , and (p"2]), respectively. 



(17) 
(18) 

(19) 
(20) 



3.2. Baryons 



In the case of cosmologies without cold dark matter, the transfer function departs from unity below 
the sound horizon to exhibit a series of declining peaks due to acoustic oscillations. The small-scale exact 
solution of equation ( |l3| ) suggests that these may be written as the product of a declining oscillatory term, 
a suppression due to the decay of potentials between the equality and drag scales, and an exponential Silk 
damping. We therefore write 



n = 



1 + (ks/5.2) 2 1 + {I3 b /ksf 



-(fe/fesiik) 1 



jo(ks). 



(21) 



Here the spherical Bessel function jo(x) = sin(x)/x is a piece that approaches unity above the sound horizon 
but oscillates below it. The envelope in square brackets traces the zero-baryon CDM case above the sound 
horizon and then breaks to a constant multiplied by an exponential Silk damping factor. We attach the Silk 
damping factor only to the second term because such diffusion can only occur on scales below the sound 
horizon, where only the second term contributes significantly. This subtlety marginally improves the fit. The 
sound horizon s, Silk scale fcsiik, and amplitude suppression were given in equations (^J), (f7|), and ( p"i| ) 
respectively; we now discuss § and /3&. 

While the nodes of the baryonic transfer function asymptotically approach those of sin(/cs) for ks ^ 1, 
the first few nodes (ks <J 10) fall at higher k than predicted by sin(ks). This shift is due to the contribution 
of the baryon density perturbation itself at the drag epoch and reflects the fact that at the sound horizon 
velocity overshoot is not the dominant effect. This effect increases with il^h 2 because the time available for 
velocity overshoot (see eq. |lj]) decreases as (^d/^cq) 1 ^ 2 and is only weakly dependent on the baryon fraction. 
We have verified this explanation of the node shift by isolating the density and velocity contributions of the 
baryons at the drag epoch from numerical evolution codes. 



We address this shifting of the nodes phenomenologically by introducing the quantity 

s 



(1 + (/?node/fcs)3) 



3^/3' 



(22) 
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For ks 3> /3 n odc, s — ► s, restoring the sinusoidal nodes. However, at ks <J /3 n odc, s rs fcs 2 //3 no de < s, moving 
the nodes to higher fc. We find 

/3„odc = 8.41(f2 ^ 2 ) - 435 , (23) 
independent of the baryon fraction. Hence the effect gets smaller at low f2o as expected. 

The amplitude a& specifies the small-scale asymptotic contribution of the velocity portion of the acoustic 
oscillation. Two effects modify this amplitude at large scales. Above the sound horizon, velocities contri- 
butions to the transfer function fall off. Furthermore, the amplitude declines if CDM dominates the energy 
density of the photon-baryon system when the wavelength enters the horizon. This occurs due the absence of 



feedback in the gravitational driving of the photon-baryon oscillator (HS96 §3.1). The net result is a cut-off 
associated with the sound horizon which moves to smaller scales as fio/i 2 increases and/or Oj, / Oo decreases. 
We describe this in equation (fell) by turning on the velocity term at the characteristic scale /?f,s, where 



Pb = 0.5 + ^ + (3 - 2^)^(17.2^)2 + i. (24 ) 

"0 "0 



3.3. Performance 

For the parameter range 0.025 <; floh 2 < 0.25 and < fib/fio < 1 the fitting formula works quite 
well. For fully baryonic models (i.e. fl c — 0) the fractional residuals are nearly always under 10%. As the 
baryon fraction decreases, the accuracy improves due to the increasing contribution of the simpler CDM 
piece. Residuals smaller than 5% are typical for f^/TJo < 0.5. Note that we quote the residuals as the 
difference between the fit and the numerical results divided by a non-oscillatory envelope that is defined by 
replacing sin(ks)/ks in equation ( ^i"| ) by [1 + (fcs) 4 ] -1 / 4 . This envelope matches the knee of the transfer 
function and grazes all the subsequent maxima. 

In Figure || we display four comparisons of the fitting formula relative to the numerical results. Also 
shown are the residuals relative to our envelope. The reason for the degradation in the fit at the shortest 
scales is that small errors in the sound horizon (c.f. the end of §|2.2| ) produce significant errors in the phase 
of the oscillations, producing order unity residuals. However, the fitting formula reproduces the correct 
amplitude and hence the Silk scale. 

The most serious systematic error in the fitting formula occurs for Q,^h 2 > 0.25. In the baryon sector of 
these cases, the drop between ks « 1 and the oscillations at ks J> 5 becomes quite precipitous. Our formula 
does not decline this quickly, causing the amplitude of the first valley to be significantly overestimated. Later 
peaks and valleys are underestimated in an attempt to compensate. One can see the beginnings of this trend 
in the flo = fi& = 1 example in Figure |c the problem gets more severe for higher fl^h 2 . 

A less important systematic effect occurs for high baryon fraction (0.7 < ilb/ilo < 0.9), low-Si models. 
Due to a small shift in the CDM break scale (eq. jl^]) that we have chosen not to model, the first valley is 
systematically underestimated by ~ 10-15% in amplitude. This problem does not extend to lower baryon 
fractions. 



4. Phenomenology 



There are a number of phenomenological trends as a function of cosmological parameters. The two 
basic effects that arise from the inclusion of baryons are the introduction of oscillations and the suppression 
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Fig. 3. — Four examples of the fit compared to numerical results. The larger plots show the numerical result 
(solid) and the fit (dashed). The smaller subplots show the residuals, defined as the difference between the 
two divided by a non-oscillatory envelope. Note that in the fully baryonic models, the oscillations have 
alternating sign in the transfer function. Also shown is the zero baryon case (dotted); note the strong 
suppression on scales below the sound horizon due to the baryons. 
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Fig. 4. — The location of the first peak in Mpc" 1 as a function of cosmological parameters. 

of power below the sound horizon, with a corresponding sharpening of the bend around the sound horizon. 
We discuss these two in turn. 



4.1. Baryon Oscillations 

Two interesting aspects of the baryonic oscillations are the location and amplitude of the peaks and 
troughs. Well under the sound horizon we expect them to fall at k = rmr/2s where m = 3, 7, 11 . . . for 



troughs and m = 5, 9, 13 . . . for peaks. However as described in §3.2, the first few oscillations are shifted 
according to the parameter /3 n odc- As flgh 2 increases, the first few peaks and troughs are progressively shifted 
to higher k. A corollary to this shift is that the ratio of the node locations becomes smaller as one raises £Iq, 
i.e. the valleys and peaks become slightly narrower. The location of the first peak is conveniently fit as 

fepeak= ^(1 + 0.217^) (25) 

where 

44.5 ln(9.83/f^ 2 ) , r 
s = v 1 ' Mpc 26 

approximates the sound horizon to ~ 2% over the range flbh 2 ;> 0.0125 and 0.025 <; floh 2 < 0.5. The value 
of fcpeak as a function of cosmological parameters is shown in Fig. ^. 

The amplitude of the oscillations also has a non-trivial dependence on the cosmological parameters. 
The oscillations of course grow stronger as the baryon fraction increases. However, at fixed baryon fraction, 



- 13 - 



(a) First Trough (b) First Peak 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 



£l h 2 Q h 2 

Fig. 5. — The fractional enhancement of power due to (a) the first valley and (b) the first peak relative to 
non-oscillatory CDM portion of the transfer function, (T/T c ) 2 — 1, at the appropriate wave vector from the 
fit. Unlabeled contours are at (a) -0.8, 0, 1, 3 and (b) 8.0, 16.0 



they are weaker compared with CDM contributions in high Slo^ 2 due to the increase in the time available 
for the CDM to grow between equality and the drag epoch. While the full series of peaks and valleys may 
be impossible to observe due to non-linear structure formation, the first valley and peak are generally in 
the linear regime. In Fig. |, we show the fractional enhancement of power due to the oscillations over the 
smooth CDM contributions. The first peak grows monotonically with the baryon fraction. The first trough 
is more subtle: while the transfer function at this location simply declines as the baryon fraction increases, 
when T(k) goes negative the power, which is the square of T, will regenerate. Perfect cancellation of the 
baryon and CDM contributions occurs along the contour labeled "-1" in Fig. ||a; above this line the trough 
in amplitude becomes a peak in power as the baryon contributions come to fully dominate. A useful rough 
scaling as to when oscillations become important is given by 

-±^Q h 2 + 0.2 (27) 

"0 

which crudely describes the region where the change in power is greater than ~ 20%. 

A related trend is the increase with £loh 2 of the sharpness of the decline in the baryonic sector from the 
knee at ks ss 1 (T « 1) into the series of oscillations below the sound horizon. This is most easily seen in 
the fully baryonic models of Figure ||; near k — 0.1ft.Mpc -1 in the VL^h 2 — 0.25 model, the transfer function 
drops a factor of ten in under half a decade in k. This break becomes even more striking in higher floh 2 
cases. 
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4.2. Effective Shape 

As we have seen, if fif,/fio <J f2 /i 2 + 0.2, the main effect of the baryons is not to introduce oscillations 
into the transfer function but to suppress the k~ 2 tail from the growth of CDM perturbations. This occurs 
both because the CDM portion T c is suppressed by a c and because the baryonic portion XJ, is providing 
essentially no power below the sound horizon. As noted just above, the latter transition can occur quite 
quickly. These suppressions indicate that the shape of the transfer function must change, with a break near 
the sound horizon. It is useful to quantify this effect. 

Let us work forward from the zero baryon case. Here the transfer function is parameterized by q cx k/k eq , 
more commonly expressed as a shape parameter T — Oq/i, where 



k 



/iMpc 



-©2. 7 /r. 



(28) 



A commonly used fitting formula to the zero-baryon limit was presented by Bardeen ct al.| (1986, eq. G3). 
However this formula fits neither the exact small-scale solution of §2.2 nor does it have the quadratic deviation 
from unity required by the theory. The latter is a fundamental requirement of causality ( [ZePdovich" 1965), 
in that one power of k must arise from stress gradients generating bulk velocity and a second from velocity 
gradients generating density perturbations. In fact the coefficient of this quadratic deviation can be calculated 
perturbatively if the stress gradients are dominated by the isotropic (pressure) term. 

The following functional form satisfies these criteria and is a better fit to the zero baryon case extrapolated 
from trace-baryon models calculated by CMBfast 



Mq) 
L (q) 
Co(9) 



Lo + C q 2 ' 
ln(2e+ 1.8. 



(29) 



14.2 



731 



l + 62.5g 



Note that this form is not only more accurate than the Bardeen et al. (1986) one but is also simpler: there are 
fewer parameters and the coefficients 1.8 and 14.2 are derived theoretically. The parameter 731, the small-g 
quadratic deviation, has been fit rather than derived to account for the small correction due to anisotropic 
stress gradients. In Fig. ^ we show a comparison of this form to numerical calculations and various fitting 
formulae in the literature. Our formula agrees with numerical calculations at the same level as different 



numerical calculations agree with each other (Hu et al. 1995), i.e. to 1% through the CMB and large-scale 
structure regimes. 



The presence of baryons has commonly been included by fitting a constant shape parameter T (Peacock 



fe Dodds 1994; Sugiyama 1995). That such an approach can work on small scales can be seen as follows. 
On small scales, the effect of the baryons is a constant suppression by the factor a c Q c /D,Q. Since the 
transfer function there is proportional to (fc/To) -2 , a rescaling of To approximates this effect. However, this 
simple rescaling of T does not properly treat the region observable through large-scale structure. Above 
the sound horizon, the baryons and CDM are indistinguishable, and so the transfer function is close to that 
parameterized by Tq = Qo n - Below the horizon, if one neglects the oscillations caused by the baryons, the 
transfer function is suppressed and roughly follows that of a rescaled T. Hence the transition around the 
sound horizon cannot be fit by a single T. 



A reasonable fit to the non-oscillatory part of the transfer function can be written by rescaling r e g(fc) 
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0.01 0.1 

k (h Mpc- 1 ) 



Fig. 6. — Comparison of various zero baryon transfer functions to the numerical calculation of CMBfast 

9j) agrees with CMBfast to 1%, whereas previous 



(Seljak & Zaldarriaga 1996). The fitting form of equation 
approximations [Pccblesj (1 982) [P82], jBtarobinsky fc Sahni| (1984) [SS84], and jBardeen et al"| (1986) [BBKS86]. 
Also shown is the form of Bond & Efstathiou (1984) [BE84], a 3% baryon fraction model that has often been 
used as a zero-baryon proxy. 



as one moves through the sound horizon 

T cS (k) = Voh(a r+ i+ \ o 2 ks)4 ), (30) 

q /n \ 2 

or = 1 - 0.328 ln(43m /i 2 ) + 0.38 ln(22.3ft /i 2 ) I -M . (31) 

Dehning q e g as in equation (p8|), we find that T{k) w Tb(q c ff). Figure ^ shows an example. Note the simpler 
form of s in equation ( |26| ) may be used here, ar is nearly a c Q c /Qo, the radicand of which is plotted in Fig. 
|2|; we provide the above form for simplicity and to account for small deviations at the higher Vl^h 2 values. 
The latter arise because the fit has been optimized in the observable region k <J 0.1 Mpc -1 , where the CDM 
transfer function is not quite in its k~ 2 small-scale limit. Of course, neglecting the oscillatory contributions 
is not a good approximation for Ob/fio 0.5. In addition to the obvious omission of the wiggles, we have fit 
ar to the transfer function T(fc) and not the observable T(k) 2 . This neglects the power introduced by the 
square of the oscillatory term. 

We caution the reader that the small-scale asymptote r c ff = ar^oh arises from the normalization 



rather than shape and therefore should not be conflated with T derived from redshift surveys (Peacock & 



Dodds 1994). Pure normalization distinctions are not observable with current redshift surveys; instead, one 
estimates T by fitting power spectra of arbitrary normalization and relying on the differences in shape and 
power-law slope introduced by the variation in T. However, when combined with the COBE normalization 



e.g. Bunn fc White) 1997), the normalized predictions from equations ( |3C| ) and ( |3l| ) are appropriate. 
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Fig. 7. — Distortions to the shape of the transfer function as measured by the true transfer function over 
the zero baryon form. The shape function of equation ( |30[ ) adequately describes the true curve apart from 
the oscillations, whereas the constant shift in T described by 
(1995) [S95] do not. 



Peacock & Dodd; 



(1994) [PD94] and jBugiyarna 



However, equation ( |30|) and Fig. fj] do make clear that a simple rescaling of T is not accurate near the 
sound horizon. Future measurements of large-scale structure expected from the Sloan Digital Sky Survey 
(Gunn & Weinberg 1995) and 2-degree Field ( Tayloij 1995) surveys will be sufficiently precise as to detect 
the deviations from a single T model (Tegmark 1997; |Goldberg fc Straus^ 1997). A key signature is that the 
bend from T sa 1 to T oc fc~ 2 becomes sharper as baryons are introduced. Detecting this change will require 
that scales around s/ir, c.f. equation (Bq), are well-observed. 



5. Discussion 

We have presented an accurate, well-motivated fitting form for the transfer function of a general CDM- 
baryon universe. The formula is generally accurate to better than 10% in fully baryonic universes and 
better than 5% in cosmologies with f^/rio ^ 0.5. While the available numerical codes will provide yet more 
accurate transfer functions, the fitting form here should be useful for characterizing trends in cosmological 
parameter space. Moreover, by separating the various physical aspects in the analytic form, we hope to 
provide physical intuition for the effects, their interrelationships, and their correspondence with analogous 
effects in the CMB. 

As applications of the form, we gave quantitative assessments of the location and amplitude of the 
baryonic features in the linear regime. We quantified the suppression of power on scales below the sound 
horizon due to the admixture of baryons and showed that this suppression is not well-fit by a rescaling of T 
if one is probing scales near the sound horizon. An alternative model, based on an interpolation in T, gives 
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a reasonable fit to the non-oscillatory portion of the transfer function provided that fifc/Tio ^o^ 2 + 0.2. 
Finally, we gave a new fitting form for the zero-baryon limit that is more accurate on small scales than those 
used previously. A summary of how to use these various formulae is given in Appendix A. 



Baryonic features, like the acoustic peaks in the CMB (Hu & White 1996), transcend the adiabatic 
CDM paradigm discussed in this paper. In fact any gravitational instability model where reionization took 
place no earlier than z c i ~ 270(£lo^ 2 ) 1 / 5 , the Compton drag epoch for a fully ionized universe, must possess 
acoustic effects at some level. Evidence of these effects is strong indication that fluctuations were generated 
in the early universe. Moreover, when combined with CMB observations, baryonic features allow strong 
consistency tests for the predicted growth of fluctuations as both sets of features reflect the underlying 
acoustic oscillations before recombination. A measurement of the sound horizon from the matter power 
spectrum combined with its angular extent from the CMB would allow an angular diameter distance test for 
curvature in the universe that is largely free of cosmological assumptions. Baryonic features in the matter 
power spectrum thus represent a valuable resource for cosmological information, but one that may be difficult 
to mine observationally. 
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A User's Guide 



In this paper, we have presented a fitting form for the transfer function in a CDM and baryon cos- 
mology with adiabatic perturbations. The formula is given as equations (^)-(^), (|l0|)— (|l2[) , and (|l4|)-(24) 



Alternatively, if one prefers a simpler form that accurately represents the baryon-induced suppression on 
intermediate scales but that ignores the acoustic oscillations in the transfer function, one may use the form 
given in equations (|2^) and (|2g|)-(|3"l|). The oscillations are fairly small (<; 20% modulation in power) for 
Qb/Oo <J Vlati 1 + 0.2, so in these cases the simpler form will be appropriate for many applications. We also 
provide in equation ( |29| ) a more accurate form for the zero-baryon transfer function. Electronic versions of 
the formulae in this paper may be found at 



http: / / www.sns.ias.edu/~whu/transfcr/transfer.htmj . 



The power spectrum of the density fluctuations is then proportional to the initial power spectrum times 
the square of the transfer function. In the most usual case, the initial power spectrum is taken to be a 
power-law, so that P(k) cx k n T 2 (k), where n = 1 is the familiar Harrison-Zel'dovich-Peebles scale-invariant 
case. 

While the transfer function is independent of late-time effects such as the presence of a cosmological 
constant or curvature, the magnitude and time-dependence of the normalization of the power spectrum 



does depend on these effects. Bunn & White (1997) calculate the present-day normalization of the power 



spectrum implied by the 4-year COBE anisotropy measurement and provide the following fitting forms for 
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flat and open cosmologies: 



jU3 / rk\ 3+n 



A 2 (fc)L =0 - ^m = S%[- m ) T\k), (Al) 
S H = 1.95 x io-5o o - - 35 -°- 191n0o - - 17 "e- fl - - 14 ™ 2 (A = 0), (A2) 

<Jff - L94x 10 -5 fi -0.785-0.051„no e -0.95fi-0.169n 2 (A = 1 - Q ), (A3) 

where h = n — 1 and the contribution of tensor perturbations to the observed anisotropies has been assumed 



to be zero (see Bunn & White 1997 for more details). The la statistical uncertainty is 7%, and the error in 



the above fits are much smaller than this for 0.2 < fl < 1 and 0.7 < n < 1.2. 

To extend this normalization to earlier times, one needs to scale the power spectrum by the square of 
the growth function D\{z) (Peebles 1980). As is well-known, for CIq = 1, D\ = (1 + z)^ 1 at redshift z. For 



other cosmologies, one can use the approximation (Lahav et al. 1991; Carroll et al. 1992) 



B x {z) = (1 + Z )^^M. \n{zfl 7 - Q A (z) + [1 + Q(z)/2][1 + fi A (z)/70]} * . (A4) 
Here Q(z) and il A {z) are the density parameters as seen by an observer at redshift z; hence 

<Xz) W + z) 3 (A5) 

nA(z) = n A + n R (i + zr + n (i + zr> (A6) 

where Qq is the matter density in units of critical, fl A = A/3Hq represents the cosmological constant A, and 
flu = 1 — fig — Q A represents the effects of curvature. 

The normalization of the power spectrum is such that the variance of mass fluctuations inside a sphere 
of radius R is 

where ji(x) = (x cos a; — sina:)/a; 2 . 
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